Perimetric complexity of binary 
digital images 

Notes on calculation and relation to visual 
complexity 

Andrew B Watson 


Perimetric complexity is a measure of the complexity of binary pictures. 
It is defined as the sum of inside and outside perimeters of the 
foreground, squared, divided by the foreground area, divided by An . 
Difficulties arise when this definition is applied to digital images 
composed of binary pixels. In this paper we identify these problems and 
propose solutions. Perimetric complexity is often used as a measure of 
visual complexity, in which case it should take into account the limited 
resolution of the visual system. We propose a measure of visual 
perimetric complexity that meets this requirement. 


■ Background 

Perimetric complexity is a measure of the complexity of binary pictures. It is defined as the sum of 
inside and outside perimeters of the foreground, squared, divided by the foreground area, divided 
by 4 7T. The concept of perimetric complexity was first introduced (and called dispersion) by At- 
tneave and Arnoult [1] in an effort to explain the apparent preceptual complexity of visual shapes. 
In the field of image processing, the concept appears as its inverse compactness [2-4]. The con- 
cept was given new life (and a new name) in 2006 by Pelli et al. , who showed that the efficiency 
of letter idetification was nearly proportional to perimetric complexity [5]. It has since become a 
popular metric in a variety of shape analysis applications, including human letter identification 
[5-7], handwriting recognition [8], evolution of graphical symbols [9], and design of graphical anti- 
spam technologies [10-12]. 

In this note we develop Mathematica functions to compute perimetric complexity of binary digital 
images, and illustrate their application. The code is comaptible with Version 8 of Mathematica. 

Although the concept of perimetric complexity is clear when applied to continuous plane shapes, 
complications arise when the concept is extended to binary digital images. We discuss these com- 
plications and suggest suitable solutions. We also introduce the concept of visual perimetric com- 
plexity, which takes into account the blurring action of the human visual system. 

We begin by illustrating the application of the function PerimetricComplexity to a binary 
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image. 


cat = 



PerimetricComplexity [cat] 

{2110.08, 89588, 3.95494} 

The output is a list containing the perimeter (in pixels), the area (in square pixels), and the complex- 
ity. In the following sections we will describe the derivation of this function, as well as the options 
that may be used to control its operation. 


■ Perimetric complexity of geometric shapes 

Perimetric complexity is a measure of the complexity of binary pictures. In a binary picture, one 
or several regions of the same color (white) are defined as foreground, and the remainder (black) 
as background. Perimetric complexity C is defined here as the sum of inside and outside perime- 
ters of the foreground P, squared, divided by the foreground area A, divided by 4 n 

P 2 

C = . (1) 

A4tc 

In the remainder of this note, unless otherwise noted, we will use the term complexity as syn- 
onomous with perimentric complexity. 

We begin with the example of a circular disk with unit radius. 

Graphics [ {White , Disk [ ] } , PlotRange -> 2 , Background -> Black] 



Here the perimeter is 2 7r, and the area is n, so the complexity is 


(2 TO 2 

C = = 1 

nAn 


( 2 ) 


It can be shown that the disk is the shape with the lowest complexity. The normalizing constant 
An in the definition leads to a unit value for this most simple shape. As a consequence, any other 
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value of complexity is easily compared to that of the circular disk. Pelli et al. [5] suggest that com- 
plexity is closely related to the number of visual features in a shape. In that sense, we could say 
that the circular disk has only a single feature. 

Our next example is a square with unit sides. 

Graphics [ {White, Rectangle [-{ 1 , l}/2, {1, l}/2]}, PlotRange 1 , 
Background -> Black] 



The perimeter here is 4, and the area is 1, so the perimetric complexity is. 

4 2 4 

C = = - « 1.27 (3) 

4n n 

If we add a square hole in the center, with sides 1/2, there will be an interior perimeter as well, as 
shown here. 

Graphics [ {White, Rectangle [-{ 1 , l}/2, {1, l}/2] 

, Black, Rectangle [- {1, l}/4, {1, 1} / 4] } 

, PlotRange -> 1, Background -> Black] 


Now the total perimeter is the sum of inner and outer perimeters, and the area is the difference in 
areas of the squares, so 

(4 + 2) 2 12 

C = = — « 3.82 (4) 

(1 - 1 /4) 4n n 

So according to this measure, the square with a hole is about three times as complex as the square. 

Some important observations about complexity are 1) it is dimensionless, 2) it is independent of 
scale or orientation, 3) it is additive. By additive we mean that the complexity of a pair of shapes, 
considered as a single shape, is equal to the sum of their PCs computed separately. 
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■ Perimetric complexity of plane curves 


Although it is beyond the scope of this note, we note for reference that if a shape is defined by a 
closed parametric curve, its exact complexity can be obtained using calculus methods [13]. Specifi- 
cally, if over an interval a < t < b the functions x(t) and y(t) and their derivatives i;(t) and y(t) are 
continuous, then the curve described has a length 


P = 



( 5 ) 


and an area 


A = - 



( 6 ) 


■ Perimetric complexity of binary digital images 

A digital image is defined here as a rectangular array of square pixels. A binary digital image con- 
tains pixel values of 1 (white) and 0 (black) only. The foreground consists of the white pixels. 

The original definition of complexity relies upon the notion of a perimeter, which has no unique 
analog in the context of digital images. However, two definitions of perimeter are available, as de- 
scribed below. 


□ Using “PerimeterLength” 

The first definition we consider is the most straightforward. Consider a binary image consisting of 
a single white pixel : 

ImagePad [Image [{ {1} }] , 1] 



It seems natural to define the perimeter of this shape as 4, and the area as 1, so C = 4 In, the same 
as the square discussed earlier. 

Now consider this shape consisting of 3 white pixels. 

ImagePad [Image [{ (1 , 0}, {1, 1}}], 1] 



Here the perimeter, consisting of the exposed pixel faces, is 8, and the area is 3, so C = 16/(3 n). 

Extending this idea, we can define the perimeter as the sum of the exposed faces of pixels in the 
foreground. 

Version 8 of Mathematica includes a set of operators from the discipline of mathematical morphol- 
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ogy. These can be used to easily calculate perimetric complexity. To illustrate this we begin with a 
binary image with several separated parts. 


yi 



The MorphologicalComponents operator finds connected regions, and labels them with inte- 
gers. The Colorize operator visualizes these regions by assigning colors to each label. The 
CornerNeighbors->False option ensures that only 4-connected neighborhoods are considered. 

MorphologicalComponents [yi, CornerNeighbors False] // 

Colorize 



The ComponentMeasurements operator returns a selected set of measurements about each re- 
gion. In this case we are interested in the area and the perimeter length. The results are returned as 
set of rules, showing the results for each region. 

cm = ComponentMeasurements [ 

MorphologicalComponents [yi, CornerNeighbors -> False] , 

{ "PerimeterLength" , "Count" } , CornerNeighbors -> False] 

{ 1 -> { 56 , 98 }, 2 -> { 198 , 469 }, 3 ^( 138 , 371 }} 

We can combine the perimeters and areas of the several regions, and then compute complexity in 
the usual way. 

Total [Last /@ cm] 

{ 392 , 938 } 

% [ [1] ] * 2 /%[ [2] ] / (4 Pi) // N 
13.0365 

The preceding calculations are implemented in the function PerimetricComplexity, defined in 
the Appendix. To obtain the previous result, we evaluate: 

PerimetricComplexity [yi , Filter -> None, 

Method -> "PerimeterLength"] // N 

{ 392 . , 938 . , 13 . 0365 } 

The option Method^'PerimeterLength" ensures that we use the definition of perimeter described 
above. The option Filter->None will be explained below. 
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For future reference, to distinguish it from variants that we will consider, we will call this the 
“raw” perimetric complexity. 


□ Using “PolygonalLength” 

The second definition of the perimeter of a connected region in a binary digital image is to con- 
sider the perimeter pixel centers as points on a lattice, and to define the length as the sum of the 
sides of the polygon defined by those points. This estimate of the perimeter is obtained from 

ComponentMeasurements by using the measurement “PolygonalLength” . 

cm = ComponentMeasurements [ 

MorphologicalComponents [yi, CornerNeighbors -> False] , 

{ "PolygonalLength" , "Count" } , CornerNeighbors -> False] 

{1 -> { 39 . 6985 , 98 }, 2 ^{ 159 . 439 , 469 }, 3 ->{ 117 . 012 , 371 }} 

Note that the measures of perimeter length are smaller than before. 

And again complexity can be easily computed: 

N [tt [ [1] ] * 2 / ti [ [ 2 ] ] / (4 Pi)] &@Total [Last /@ cm] 

8.47953 

This variant of perimetric complexity is implemeted with the following options: 

PerimetricComplexity [yi , Filter -> None, 

Method -> "PolygonalLength"] / / N 

{ 316 . 149 , 938 . , 8 . 47953 } 

Or, since Methods ’’PolygonalLength” is the default, 

PerimetricComplexity [yi , Filter -> None] // N 

{ 316 . 149 , 938 . , 8 . 47953 } 


■ Approximating complexity of continuous shapes 

We introduced the concept of perimeric complexity with a few continuous shapes, such as a 
square and a circle. In these cases, complexity is easily calculated, because we have simple formu- 
las for the area and perimeter. It might be imagined that complexity of the continuous shape could 
be approximated by computing the complexity of a discrete sampled image, rendered from the 
shape. As we shall see, this assumption is not strictly correct. 

Consider the circular disk. As noted at the beginning, it has C = 1. 
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disk = Graphics [ {White, Disk[]}, Background -+ Black] 



Now we consider an image rendered from the continuous shape. We render it into a certain size im- 
age. We change the color so that the foreground is white, as is our convention. 

size = 2*7 + 1; 

diskimage = Image [ disk, "Bit", ImageSize -+ {1, 1} size, 

ColorSpace "Grayscale"] 



If we compute the complexity, we find that it is too large, by 62%, relative to the continuous shape. 

PerimetricComplexity [diskimage. Filter -+ None, 

Method -+ "PerimeterLength" ] / / N 

{ 492 . , 11 873 . , 1 . 62241 } 

The reason is that the sampled image actually is more complex than the continuous shape. Its con- 
tour is jagged, while that of the continuous shape is smooth. It might be imagined that this could 
be remedied by increasing the resolution of the rendering. Here we show that belief is misplaced. 
We render at several sizes, and plot the results. Size has little effect, and the complexity never ap- 
proaches the value of 1 corresponding to the continuous shape. 
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points = Table [ 

{ size = 2 * k + 1 , 

PerimetricComplexity [ 

Image[disk, "Bit", ImageSize ->• { 1 , 1} size, 

ColorSpace "Grayscale"] 

, Filter ->■ None, Method "PerimeterLength" ] [ [3] ] } 

, {k, 4, 12}]; 

ListLogLinearPlot [points, PlotRange -> {1, Automatic}, 

Frame -> True, FrameLabel -> {"Size" , "Complexity"} , Joined -> True] 



This problem can be somewhat ameliorated by using the PolygonalLength measure of perime- 
ter length. Rather than the pure approach of measuring the exposed face of each foreground pixel, 
this measures the length of the contour that travels between the centers of those pixels. 

PerimetricComplexity [diskimage. Filter -> None] 

{ 403 . 647 , 11 873 , 1 . 09203 } 

Now the difference is reduced to 9%. Here again, the reader might think that this difference could 
be reduced to zero by enlarging the resolution (number of pixels) in the rendered image, but this is 
not so. We leave that as an exercise for the reader. The error can never be zero, because the path be- 
tween pixel centers must always be vertical, horizontal, or diagonal, so it can never smoothly fol- 
low the true circular contour. Put another way, it has a higher fractal dimension than the circle, 
and thus greater length. 


■ Pelli algorithm 

Pelli et al. [2] proposed a method for computing complexity which we quote here in full: 

The ink area is the number of l’s. To measure the perimeter we first replace the image 
by its outline. (We OR the image with translations of the original, shifted by one pixel 
left; left and up; up; up and right; right; right and down; down; and down and left; and 
then bit clear with the original image. This leaves a one-pixel-thick outline.) It might 
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seem enough to just count the l’s in this outline image, but the resulting “lengths” are 
non- Euclidean: diagonal lines have “lengths” equal to that of their base plus height. In- 
stead we first thicken the outline. (We OR the outline image with translations of the origi- 
nal outline, shifted by one pixel left; up; right; and down.) This leaves a three-pixel- 
thick outline. We then count the number of l’s and divide by 3. 

This method can be implemented using the Dilation operator, as we show here. With verbose 
=True, it shows two images: the perimeter in red, and the thickened perimeter. It returns the 
length of the perimeter, the area, and the complexity. 

2/28/11 12:10:52 (Local) ln[2]:= 

PelliMethod [image_, verbose_: False] := Module [ 

{tmpO, tmpl, tmp2, tmp3 , perimeter , area}, 

tmpO = ImagePad [ image , 2 ] ; 

tmpl = Dilation [tmpO, BoxMatrix[l] ] ; 

tmp2 = ImageSubtract [tmpl, tmpO] ; 

tmp3 = Dilation [tmp2 , CrossMatrix [1] ] ; 

If [verbose, ( Print @ GraphicsRow [ 

(ColorCombine [ {tmpl, tmpO, tmpO} ] , tmp3} ] ) ] ; 
perimeter = Total [ImageData [tmp3] , 2] /3; 
area = Total [ImageData [tmpO] , 2]; 

{perimeter, area, perimeter^2 / area / (4 Pi) } 

] 

Applying this to the three-component Chinese character, we get: 

PelliMethod [yi. True] // N 



{ 326 . 333 , 938 . , 9 . 03463 } 

The exact method gives the following: 

PerimetricComplexity [yi , Filter -> None, 

Method ->• "PerimeterLength" ] / / N 

{ 392 . , 938 . , 13 . 0365 } 

We see that the perimeter is substantially underestimated by the Pelli method in this case. This 
method has other limitations. It effectively assumes regions that are large in pixel dimensions. For 
example, consider the case of a single pixel object. As noted above, it has C = 4/n « 1.273. But the 
Pelli method yields an complexity value more than three times too large. 
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test = ImagePad [Image [{ {1} } , "Bit"], 3] 



PelliMethod [test. True] // N 



{ 7 . , 1 . , 3 . 8993 } 


■ Visual perimetric complexity 

Much of the motivation for the use of perimetric complexity is the hope that it might provide an ap- 
proximate measure of the visually perceived complexity of shapes. But this only makes sense if 
the shape is actually visible. Consider the difference between the continuous circular disk, and its 
sampled image, as discussed above. They have different perimetric complexities, no matter how 
high the resolution of the sampled version. But of course, at a certain viewing distance, they will 
be indistinguishable. 


□ Filtering 

Here we propose an approach to dealing with this problem. The idea is to first blur the image, in a 
manner consistent with visual blur, and then compute perimetric complexity. We begin with the ex- 
ample of a Chinese character. 


jun 



To make things simple, we will use Gaussian blur, although this is not an accurate description of 
human visual blur. Later we will show a more accurate form of blur. 

First we pad the image, so that the blur is contained. 
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tmpO = ImagePad [ jun, 5] 



Next we magnify the image. This allows greater flexibility in the filtering. 

mag = 4; 

tmpl = ImageResize [tmpO , mag ImageDimens ions [tmpO] ] 



Then we blur the image, in this case by a Gaussian filter with a radius of 8 pixels. 
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tmp2 = GaussianFilter [tmpl , 8] 



Then we binarize the image. Unfortunately, this requires some method of setting the threshold. 
Here we use a fixed threshold of 0.5. 

tmp3 = Binarize [tmp2 , 0.5] 



And because we imagine that the image is viewed at such a distance that the pixels are not re- 
solved, we use the PolygonalLength method. 
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PerimetricComplexity [tmp3 , Filter -> None, 
Method -> "PolygonalLength" ] // N 

{ 1944 . 82 , 24 044 . , 12 . 5182 } 

We can compare this to the unfiltered “raw” version: 

PerimetricComplexity [jun. Filter -> None, 

Method -> "Per imeterLength" ] // N 

{ 606 . , 1467 . , 19 . 9207 } 

The filtered version has substantially lower complexity, as we expect. 


□ Visual filtering using a Gaussian 

For the filtering to approximate visual blur, it must be based on the size of the original image and 
its distance from the viewer. Obviously, as the shape becomes smaller or farther from the ob- 
server, its details are more blurred, less visible, and contribute less to the visual complexity. 

The challenge is to determine the appropriate value of the Gaussian filter radius for a given view- 
ing distance. From measurements of visual sensitivity, we know that visual Gaussian blur has a 
standard deviation of about B = 0.01549 degrees of visual angle [11]. But we need to convert this 
into a radius in pixels. Recall that the image may be magnified by M before filtering. If we express 
the viewing distance V in terms of pixels (before magnification), then the size of each pixel S 
(after magnification) is approximately 

57.3 

5 = deg (7) 

MV 

The constant 57.3 is an approximation to l/ArcTan[l°]. 

1. / ArcTan[l °] 

57.3016 

By default, the radius is twice the standard deviation. So the filter radius should be 
2 B 2 BMV 

R - — = pixel (8) 

5 57.3 

Consider an example. Suppose that charimage is displayed on a typical computer screen with a res- 
olution of 72 pixels/inch, and viewed at a distance of 12 inches. Then 

V = 12 x 72 

864 

And recall that 

B = 0.01549; 
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Suppose further that the magnification is 

M = 4; 


Then 


R = 2BMV / 57.3 

1.86853 

proceding now with the steps outlined above for filtering, 

tmpO = ImagePad [ jun, 1] ; 

tmpl = ImageResize [ tmpO , M ImageDimensions [tmpO] ] ; 
tmp2 = GaussianFilter [tmpl, R] ; 
tmp3 = Binarize [tmp2 , 0.5] 

PerimetricComplexity [tmp3 , Filter -> None] // N 



{2365.71, 23 472. , 18.9741} 

Notice that the final result is not that different from the raw value. 

PerimetricComplexity [jun. Filter -> None, 

Method -> "PerimeterLength" ] // N 

{606. , 1467. , 19.9207} 

This is because the image is viewed close enough that the blur has little effect. But if we repeat the 
process, with 4 times the viewing distance, we will see that the complexity is substantially reduced. 

1 / ArcTan [1 . " °] 

57.3016 
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2.33 / 60 

0.0388333 

V = 4 x 12 x 72; 

M = 4; 

R = 2BMV / 57.3 
tmpO = ImagePad [ jun, 1] ; 

tmpl = ImageResize [ tmpO , M ImageDimensions [tmpO] ] ; 
tmp2 = GaussianFilter [tmpl, R] ; 
tmp3 = Binarize [tmp2 , 0.5] 

PerimetricComplexity [tmp3 , Filter -> None] // N 

7.47413 



{ 1987 . 85 , 23 932 . , 13 . 1395 } 

The preceding steps of padding, magnification and filtering are built in to the function 
PerimetricComplexity, as shown below. With Verbose->True, the function also shows the 
original, the filtered version, the binarized version, and the original (in red) with the perimeter of 
the filtered version (in white - within the original- and aqua - outside the original). 
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PerimetricComplexity [ jun, Magnification ->• 4 , Filter "Gaussian" , 
ViewingDistance -> 4 x 12 x 72, Verbose -> True] // N 



{ 1987 . 85 , 23 932 . , 13 . 1395 } 

The Gaussian is parameterized by a scale in degrees of visual angle. The default value is 2.33/60 
degrees. The user can experiment with different values via the option GausianScale. 


□ Accurate visual filtering with Sech 

Visual blur is more accurately represented with filters other than a Gaussian. In one simple form, 
the kernel is a Sech (Hyperbolic Secant) function [14,15]. This filter can be selected with an op- 
tion (the default) in the function PerimetricComplexity: 

PerimetricComplexity [jun. Filter "Sech" ] // N 

{ 1016 . 14 , 5977 . , 13 . 7472 } 

Or, because it is the default, 

PerimetricComplexity [jun] // N 

{ 1016 . 14 , 5977 . , 13 . 7472 } 

The Sech is parameterized by a scale in degrees of visual angle. The default value is 2.16/60 de- 
grees [14, 15]. The user can experiment with different values via the option SechScale. 

To the advanced reader, it may be of interest that the Sech kernel corresponds to a Sech Contrast 
Sensitivity Function, with a scale that is the inverse of the kernel. For example the CSF correspond- 
ing to the default scale is plotted here: 
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LogLogPlot [Sech [Pi u 2 . 16 / 60] , {u, 1, 60}, PlotRange { . 0099 , 1.01}, 
FrameLabel -> { "Spatial frequency (cycles/degree) " , 

"Contrast sensitivity"}] 



□ Accurate visual filtering with an arbitrary point spread function 

In real human eyes, blur results from both low-order aberrations such as defocus and astigmatism, 
and also from higher-order aberrations. In this example, we use a blur function defined by an array 
of values representing the filter kernel. This example is an actual estimate of the point-spread func- 
tion for an individual human observer, as measured using a device called a wavefront aberrometer, 
that includes both low and high order aberrations [ 16 ]. 

ImageAdjust [Image [psf ] ] 



We can use this blur kernel by supplying it directly to the Filter option. 
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Per imetricComplexity [char image. Verbose True, Filter -»psf] // 
N 



{ 120 . 569 , 950 . , 1 . 21768 } 

We will not adress this topic in detail, but when a kernel is supplied, the pixels of the kernel and of 
the magnified image must be of the same size, in degrees of visual angle, for the filter to be accu- 
rate. The size of the magnified pixels in degrees will depend upon the viewing distance and the 
magnification, as described above. Here we magnify the image, and get a different result. 


Per imetricComplexity [char image. Verbose -> True, Filter -> psf , 
Magnification -> 4] // N 



{ 1405 . 26 , 25 663 . , 6 . 12344 } 


□ Magnification parameter 

The Magnification parameter should be set with a value that ensures that a filter kernel, if pre- 
sent, has enough pixels in it to adequately represent the filter. For the Sech filter, we generally 
want a width of least 3 times the scale, and at least 8 pixels. This means that the magnification 
should be greater than 

8x57.3 

M > (9) 

3 V B 

This rule is effectively implemented by the default Magnif ication->Automatic option. 


□ Binarization 

After applying visual blur, it is necessary to binarize the image before calculation of perimetric 
complexity. There are many ways to binarize an image, and mathematica offers many of them as 
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options. The simplest is to use a fixed threshold. Since our images are initially defined as 0 or 1, a 
natural choice of threshold is 0.5. One drawback of this choice is that as images become severely 
blurred, no pixels may remain that exceed the threshold. From a perceptual point of view, the 
mean might appear a reasonable choice. As the image blurs, all pixels revert towards the mean, but 
some always remain above the mean until a uniform image is reached. A drawback of the mean is 
that it is influenced by the area of the background. For this reason, we adopt the fixed value of 0.5 
as the default threshold for binarization. 

We should acknowledge that a more valid visual thresholding scheme might be devised, that better 
reflects our perceptual segregation of areas into light and dark. This is a topic for future research. 

It should also be acknowledged that for severely blurred images, considerable grayscale informa- 
tion remains that is lost in binarization. Thus we should question whether perimetric complexity is 
an appropriate measure for such images. 

To illustrate the problem, we show an example of a character viewed at 10 feet on a display with 
100 pixels/inch. Note that the blurred image displays internal grayscale structure that is not con- 
veyed by the binarized version. 


Per imetricComplexity [char image, ViewingDistance -> 10 x 12 x 100, 
Verbose ->■ True] // N 



{ 199 . 51 , 1604 . , 1 . 97476 } 


□ Viewing distance 

For a given shape, complexity will decline with viewing distance, as a result of visual blur. Here 
we illustrate this effect with an example Chinese character. First we find the raw complexity, 

raw = 

PerimetricComplexity [ jun. Filter -> None, 

Method -> "PerimeterLength" ] [ [3] ] // N 

19.9207 


Now we compute complexity for viewing distances ranging from 3 inches to 10 feet, assuming a 
display with a resolution of 100 pixels/inch. For reference, we show as red lines the raw complex- 
ity, and the theoretical limit of 1 (a circular disk). As it should, the visual complexity proceeds 
from one of these limits to the other as distance increases. 
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{min, max} = {3, 10x12}; 
points = 

{tt, PerimetricComplexity [ jun, ViewingDistance -> tt 100] [[3]]} & /@ 
Range [min, max, 3] // N; 

ListPlot [points , Frame ^ True, Joined True, 

Epilog -> {Red, Line [{{min, raw}, {max, raw}}]. 

Line [{{min, 1}, {max, 1}}]}, 

PlotRange ->• {{0, Automatic}, {0, 22}} 

, PlotStyle ->■ Blue, Axes -> False, 

FrameLabel -> { "Distance (inches)", "complexity"}] 



At very small viewing distances (in pixels) the blur has little effect on each pixel, so the visual com- 
plexity approaches the raw value. As a rule of thumb, this asymptote is approached when the size 
of each pixel exceeds 1/4 degree. 

It is reassuring to know that the algorithm does approach the correct asymptote as distance in- 
creases. Here we show the intermediate images for a case of extreme blur (distance = 30 feet). 


PerimetricComplexity [char image. Verbose True, 
ViewingDistance -4 30 x 12 x 100] 



{ 104 . 669 , 818 , 1 . 06579 } 

But a note of caution is warranted. Consider the effect of distance on our other example character 
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). We show it here along with the previous plot. Note that the curves cross, so that at large 


distances (large blurs) the “simpler” character becomes the more complex of the two. 


points2 = 

{tt, PerimetricComplexity [yi , ViewingDistance -> tt 100] [[3]]} & /@ 
Range [min, max, 3] // N; 

ListPlot [ (points, points2}. Frame ^ True, Joined -> True, 

Epilog -> (Red, Line [{{min, raw}, {max, raw}}]. 

Line [{{min, 1}, {max, 1}}]}, 

PlotRange -> {{0, Automatic}, {0, 22}} 

, PlotStyle {Blue, Green}, Axes False, 

FrameLabel -> { "Distance (inches)", "complexity"}] 



Q i i i I i i i I i i i I i i i I i i i I i i i LJ 

0 20 40 60 80 100 120 


Distance (inches) 


This makes sense, since the densely packed features of the “complex” character blur onto each 
other , while the more widely separated features of the “simple” character remain distinct. This is il- 
lustrated in the following, which shows the intermediate images for the two characters when 
highly blurred (distance = 8 feet). 
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PerimetricComplexity [tt, ViewingDistance 8 x 12 x 100, 
Verbose True] & /@ {charimage, charimage2}; 



But the conclusion we must draw is that even the relative complexity of different shapes cannot be 
known without specifying the viewing distance. 

As a default value, we use the quantity ViewingDistance->48 /ArcTan[ 1°] . This corresonds 
to a visual resolution of 48 pixels/degree. It is commonly encountered resolution, about that 
achieved by a display with 100 pixels/inch viewed from 27 inches. But in actual use, it is advised 
to use the actual viewing distance rather than relying on this default. 

Here is a panel that allows you to experiment with different viewing distances. The complexity, 
and the diagnostic images are shown. 

Manipulate [ 

{«[ [3] ] // N, GraphicsRow [tt [ [4] ] , ImageSize 400] } &@ 
PerimetricComplexity [charimage, Returnlmages -4 True, 
ViewingDistance inches 100] 

, {inches, 1, 200}] 


■ Recommended practice 

The reader, and potential user of PerimetricComplexity may be daunted by the number of vari- 
ants of the operator, or unsure what parameters to use. Here we offer suggestions that simplify this 
matter. 

In general, we recommend using the function with default parameters. As a reminder, these are: 
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Options [PerimetricComplexity] 

{Magnification Automatic, Pad -> Automatic, 

Verbose -> False, Threshold ^ 0 . 5 , Method PolygonalLength, 

Filter ^ Sech, Normalized True , ViewingDistance -» 2750 . 48 , 
MorphologicalOperators True, SechScale -> 0.036, 

GaussianScale -> 0.0388333, Returnlmages -> False} 

The only option that should be specified for the typical use of this function is ViewingDistance. 
This specifies the distance from the eye to the image, specified in uinits of pixels. The default is 
ViewingDistance->2750.48, consistent with a display having 96 pixels/inch, viewed at 28.65 
inches (equal to a display with 48 pixels/degree of visual angle). 

It is difficult to imagine a case in which vision, in some form, would not be used to view the shape 
in question. If such a case arises, however, the “raw” complexity can be measured with the follow- 
ing options: 

SetOptions [PerimetricComplexity, Filter -> None, 

Method -> "PerimeterLength" ] 

This will also be an appropriate measure when the pixels are very large (larger than 1/4 degree). 

When trying to approximate the raw complexity of a continuous shape by means of a sampled rep- 
resentation (e.g., a circle via an image of a circle), the following options will yield the lowest er- 
ror. But as noted above, the error will still be significant. 

SetOptions [PerimetricComplexity, Filter -> None] 


■ Examples 

We conclude with two examples of the application of PerimetricComplexity . 

The first example is a set of three binary images. Below each image we print the complexity. 

pictures = (cat, horse, family}; 
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Grid [ 

{pictures. 

Style [NumberForm [Last [PerimetricComplexity [tt] ] , 3], 18] & /@ 
pictures} 

, ItemSize 10, Spacings 2] 



3.95 8.62 22.1 


The second example is an array of characters. This array was created as part of an experiment on 
the effect of complexity on visual acuity [6]. The first row are the Sloan letters, a well known set 
of letter acuity targets [5]. The remaining six rows are sets of Chinese characters, selected so as to 
be of equal complexity within a row, but increasing in complexity from row to row [6]. The metric 
of complexity used for selection was different from that developed in this paper. 

GraphicsGrid [chararray] 


cTd]H|K|N|Q|R|S|V|Z 


VA lLi \'P + * A ^ 


n £ \Rmm ei isw ff 




m m m * w a * * 




9&m\mw 


We first apply the function to each character, and plot the results, with a different curve for each 
set. 
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2/28/11 12:10:54 (Local) In [26]:= 

c = Map [Last [PerimetricComplexity [tt] ] & , chararray, {2}]; 
ListLinePlot [c, PlotRange { { . 7 , 10.3}, {0, 18}} 

, FrameLabel -> {None, "Complexity"}] 

2/28/11 12:10:55 (Local) Out[27]= 



The figure shows that there is considerable variation in each set. If we take the mean of each set 
we see a progression in complexity, except for the last three sets. In this figure we show at the bot- 
tom an exemplar from each set. 

2/28/11 12:11:08 (Local) In [29]:= 

ListLinePlot [Mean /@ c. Mesh -> Full, PlotRange -> {{. 7 , 7.3}, {0, 15}} 


Epilog 

{Table [Inset [chararray [ [k, 1]], {k , .1}, Scaled[{.5, 0}], .6], 
{k, 7}]} 

, FrameLabel {None, "Complexity"}] 


2/28/1 1 12:1 1:08 (Local) Out[29]= 

14 

12 

^ 10 


x 

a 8 
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■ Conclusions 

We have illustrated several different methods for computing the perimetric complexity of binary 
digital images. These methods differ in how they compute the perimeter, and in whether the image 
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is blurred and binarized before the complexity calculation. We have introduced the concept of vi- 
sual perimetric complexity, and argued that in general it requires blur for a sensible estimate. We 
have described several methods of implementing visual blur. 

The computed value of visual perimetric complexity depends somewhat upon details of the calcula- 
tion, such as the presumed magnitude and nature of visual blur, and the binarization threshod. In 
this regard, we have proposed a set of standard default settings and procedures for calculation of vi- 
sual perimetric complexity. 

We have also made the observation that visual perimetric complexity cannot be estimated without 
specifying the resolution of the display and the viewing distance. As a general rule, the visual peri- 
metric complexity will approach the raw complexity when the width of a pixel exceeds 1/4 degree 
of visual angle. 


■ Appendix 

□ Functions 

Here we define several functions based on the derivations presented above. 

m PerimetricComplexity 

Clear [PerimetricComplexity] 
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PerimetricComplexity : : usage = 

"PerimetricComplexity [image_,opts Rule] Compute the 

perimetric complexity of a binary image. Perimetric 
complexity is defined as the square of the sum of 
inner and outer perimeters divided by the foreground 
area, divided by By default, foreground pixels 

are white (1) . The function returns a list: 

{perimeter, area, complexity}. Optionally, the 
image can be filtered and thresholded before 
calculation of complexity, to provide a a better 
estimate of the visual complexity of the image. 

Options available are Pad (width of padding to add 
around image) , Magnification (factor by which to 
magnify each pixel before calulation of complexity) , 
Method (\"PerimeterLength\" or \ "PolygonalLength\ " , 
the method used to compute the perimeter) , Threshold 
(a numerical threshold, or the method to use in 
thresholding operations after filtering) , and Filter 
(None: no filter, \ "GaussianX " for Gaussian filter; 
V'SechV' for a Hyperbolic Secant filter; or an 
array supplied directly as the filter kernel. For 
the gaussian and the Sech, the size of the filter 
is determined by the ViewingDistance option) , 
ViewingDistance (in pixels, used to determine the 
radius of Gaussian and Sech filters) , 
MorphologicalOperators (whether to use Mathematica ' s 
built in operators) , Normalized (whether to normalize 
complexity by the value for a disk, 47 r) . Defaults 
are: { {Magnif ication-*Automatic, Pad-*Automatic, 
Verbose-*False, Threshold-^ . 5 , Method-A "PolygonalLength\ " , 
Filter-A " Sech\ " , Normal ized-*True , 
ViewingDistance-*48/ArcTan [1°] , 

MorphologicalOperators-*True} , SechScale-* 2 . 16/60, 
GaussianScale-*2 . 33/60, ReturnImages-*False} " ; 

Options [PerimetricComplexity] = 

{Magnification -* Automatic, Pad -* Automatic, Verbose -* False, 
Threshold -* 0 . 5 , Method-* "PolygonalLength" , Filter-* "Sech", 
Normalized -* True, ViewingDistance -* 48 / ArcTan[l. °] , 
MorphologicalOperators -* True, SechScale -*2.16/60, 
GaussianScale -* 2.33 / 60, Returnlmages -* False} ; 

2/28/11 12:10:52 (Local) ln[5]:= 

PerimetricComplexity [image_, opts Rule] := Module [ 

{tmp, tmpO, tmpl, tmp2 , tmp3 , tmp4, tmp5, plength, area, 
mag, radius, pad, verbose, fdim, pc, norm, threshold, 
method, filter, vdist, pixelsperdegree, sechpixels, 
sechdegrees, morpho, sechscale, gscale, anglescale, 
r image s , image s } , 

{mag, pad, verbose, threshold, method, filter, norm, vdist. 
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morpho, sechscale, gscale, rimages} = 

{Magnification, Pad, Verbose, Threshold, Method, Filter, 
Normalized, ViewingDistance, MorphologicalOperators , 
SechScale, GaussianScale, Returnlmages} /. {opts} /. 
Options [PerimetricComplexity] ; 
anglescale = ArcTan[l. Degree]; 

If [mag = = = Automatic , mag = Switch [ filter 
, None, 1 
, _?ListQ, 1 

, "Gaussian", Ceiling [2. / (vdist anglescale gscale) ] 

, "Sech" , Ceiling [3. / (vdist anglescale sechscale) ] 

]]? 

pixelsperdegree = mag vdist anglescale; 
sechpixels = Ceiling [pixelsperdegree 3 sechscale] ; 

If [pad = = = Automatic, 
pad = If [filter = = = None, 2, 

Ceiling [pixelsperdegree sechscale] ] ] ; 

If [EvenQ [sechpixels] , sechpixels + + ] ; 
sechdegrees = sechpixels / pixelsperdegree; 
radius = 2 gscale pixelsperdegree / Sqrt [2 Pi] ; 
tmpO = Binarize [image] ; 

tmpl = If [mag > 1, ImageResize [tmpO, mag ImageDimensions [tmpO] ] , 
tmpO] ; 

tmp2 = Image Pad [tmpl , pad]; 
tmp4 = Switch [filter 
, None , tmp2 

, _?ListQ, ImageConvolve [tmp2 , filter] 

, "Gaussian", GaussianFilter [tmp2 , radius] 

, "Sech", 

ImageConvolve [tmp2 , SechKernel2D [sechpixels { 1 , 1}, 
sechdegrees {1, 1}, sechscale]] 

]? 


tmp5 = If [ImageType [tmp4] = = = "Bit", tmp4. 

If [NumericQ [threshold] , Binarize [tmp4 , threshold]. 

Binarize [tmp4 , Method -> threshold] ] ] ; 

{plength, area} = If [morpho. 

Total [ 

Last /@ 

(tst = ComponentMeasurements [MorphologicalComponents [tmp5] , 
{method, "Count" } , CornerNeighbors False] ) ] 

t 

{PerimeterLength [tmp5 , Switch [method, "PerimeterLength" , 

PixelBorderLength, "PolygonalLength" , PixelPathLength] ] , 
Total [ ImageData [ tmp5 ] , 2 ] } ] ; 

(*Print [ "mag = " ,mag, " ppd = ", pixelsperdegree, 

" sechpixels = ", sechpixels , " pad = ", pad];*) 
tmp6 = MorphologicalPerimeter [tmp5 , CornerNeighbors False] ; 
images = {tmp2, tmp4, tmp5, ColorCombine [ {tmp2 , tmp6, tmp6}]}; 
If [verbose, ( 
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Print @ GraphicsRow [images , ImageSize -> 400] ; 

)]? 

{plength, area, (plength * 2 / area) / If [norm, 4 Pi, 1]} // 
Iffrimages, Append[tt, images], tt] & 

] 

Here is an example. 

Per imetricComplexity [char image. Verbose -> True] // N 



{ 1016 . 14 , 5977 ., 13 . 7472 } 


m PerimeterLength 

2/28/11 12:10:52 (Local) ln[6]:= 

PerimeterLength: : usage = 

" Per imeterPathLength [image_, verbose_: False] Given a 
binary image that represents a set of perimeters, 
defined by connected white pixels, compute the 
total length of the perimeters. It first locates 
all foreground pixels. It then extracts the 3x3 
neighborhood of each, and using PixelPathLength, it 
computes half the distance to each of the pixels 
two nearest foreground neighbors. The perimeter is 
the sum of all these distances. If verbose is true, 
it shows a tally of the neighborhoods and their 
corresponding path lengths . " ; 
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2/28/11 12:10:52 (Local) ln[7]:= 

PerimeterLength [image_, method_: PixelBorderLength, 
verbose_: False] := 

Module [ {positions, neighborhoods, tally, perim} , 
perim = Perimeter [image] ; 

positions = Position [ImageData [perim] , 1] ; 
neighborhoods = 

ImageTake [Switch [method, PixelPathLength, perim, 
PixelBorderLength, image] , 

Sequence @@ Transpose [ { it - 1 , tt+1}]] & /@ positions; 
tally = Transpose [Tally [neighborhoods] ] ; 

If [verbose, ( 

Print ["pixels = ", Length [positions] ] ; 

Print @ 

TableForm [Transpose [ {tally [ [ 1] ] , tally [ [2] ] , 
method /@ tally [ [1] ] } ] ] ? 

)]? 

Total [tally [ [2] ] (method /@ tally[[l]])] 

] 

Here is an example. 


test = ImagePad [Image [{ {0, 1}, {1, 1}}, "Bit"], 1] 
PerimeterLength [test , PixelBorderLength, True] 



u PixelPathLength 

PixelPathLength: : usage = 

"PixelPathLength [image_] Given a 3x3 binary image, finds 
the coordinates of white pixels not at the center, 
and returns their mean distance from the center."; 
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2/28/1 1 12:10:53 (Local) ln[9]:= 

PixelPathLength [image_] := Module [ {tmpl , tmp2 , tmp3}, 

tmpl = DeleteCases [Position [ImageData [image ] , 1] , { 2 , 2}]; 
tmp2 = (tt- {2 , 2}) & /@ tmpl; 
tmp3 = Sort [Norm /@ tmp2] ; 

Mean [Take [tmp3 , 2]] 

] 

Here is an example. 


• PixelBorderLength 

2/28/11 12:10:53 (Local) ln[10]:= 

PixelBorderLength: : usage = 

"PixelBorderLength [image_] Given a 3x3 binary image, 
counts the number of black neighbor (4-connected) 
pixels. This is a measure of the length of the 
exposed border of the foreground (white) pixel."; 

2/28/11 12:10:53 (Local) ln[11]:= 

PixelBorderLength [ image_] : = 

Total [ 1 - Flatten [ImageData [image] ][[ {2 , 4, 6, 8}]]] 

Here is an example. 


Framed [t st = Image [{ {0, 0, 0}, {0, 1 , 1}, {1, 1 , 1}}, "Bit"]] 
PixelBorderLength [tst] 


tst = Image [{{0, 0, 0}, {0, 1, 1}, {1, 0, 0}}, "Bit"] 
PixelPathLength [tst] 





2 
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• SechKernel2D 

SechKernel2D: : usage = 

"SechKernel2D [samples_,degrees_, scale_] Compute a 
convolution kernel defined by a Hyperbolic Secant 
(Sech) function of distance from the origin, samples 
is a list {height , width} of dimensions of the kernel 
in pixels, and degrees is a list of the corresponding 
dimensions in degrees of visual angle, scale defines 
the width of the kernel in degrees. The origin is 
defined as Floor [sampIes/2 ] . The kernel is normalized."; 

2/28/11 12:26:19 (Local) ln[38]:= 

SechKernel2D [samples_, degrees_, scale_] := 

(tt / Total [tt, 2] ) &@ 

Sech [ 

Pi Array [N [Norm [ {tttt} degrees / samples] ] &, samples, 

-Floor [samples / 2] ] /scale] 

Here is an example. 

ImageAdjust [Image [SechKernel2D [ {32 , 32}, {1, 1}, -4]]] 



u Perimeter 

2/28/11 12:10:53 (Local) ln[13]:= 

Perimeter [image_] := Module [ {padded, positions}, 
padded = ImagePad [image, 1] ; 

positions = Select [Position [ ImageData [padded] , 1], 
PixelBorderLength [ ImageTake [padded , 

Sequence @@ Transpose [{tt-1, tt + 1} ] ] ] >0&]; 

ImagePad [ 

Image [ReplacePart [ImageData [padded] 0, tt -> 1 & /@ positions] , 
"Bit"], -1]] 


□ Initializations 

2/28/11 12:10:53 (Local) ln[14]:= 

SetOptions [Graphics , ImageSize -> 128] ; 

2/28/11 12:10:53 (Local) ln[15]:= 

SetOptions [ListLinePlot , Axes -> False, 

BaseStyle -> {10, FontFamily -> "Helvetica", AbsolutePointSize [4] } , 
Frame -> True, ImageSize -> 300] ; 
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2/28/11 12:10:53 (Local) ln[16]:= 

SetOptions [LogLogPlot , AspectRatio ->■ Automatic, Axes -> False, 
BaseStyle {10, FontFamily -> "Helvetica"}, Frame -> True, 

ImageSize 300] ; 

2/28/11 12:10:53 (Local) ln[17]:= 

SetOptions [ListLogLinearPlot , AspectRatio ->■ Automatic, 

Axes False, BaseStyle {10, FontFamily "Helvetica"}, 

Frame -> True, ImageSize -> 300] ; 

2/28/11 12:10:53 (Local) ln[18]:= 

SetOptions [ListLogLinearPlot , AspectRatio -> 1, Axes -> False, 
BaseStyle {10, FontFamily -> "Helvetica"}, Frame -> True, 

ImageSize -> 300] ; 

The next cell, containing the definition of the pointspread function, is closed to suppress printing/ 


2/28/11 12:10:53 (Local) ln[20]:= 



2/28/11 12:10:53 (Local) ln[21]:= 


yi = 


ia 


2/28/11 12:10:53 (Local) ln[22]:= 

chararray = 


o r s|v z 
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2/28/11 12:10:54 (Local) ln[24]:= 


horse 



2/28/11 12:10:54 (Local) In [25]:= 


family = 



□ Computing Complexity without using Mathematica’s Morphological 
Operators 

Some readers may wish to write code in other languages to compute complexity. For that reason 
we provide an explanation here of how to compute complexity without using Mathematica’s mor- 
phological operators. This amounts to finding aletrnate methods for computig the perimeter. These 
are incorporated in the function PerimeterLength, defined above, and derived below. These 
functions can be exercised from within PerimetricComplexity by selecting the option 
MorphologicalOperators->False. This is useful mainly for testing. 

Consider the following binary image. 
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jun 



The foreground area is easily obtained, since it is just the sum of all the white pixels: 

Total [ImageData [jun ] , 2] 

1467 


m Using “PerimeterLength” 

As noted above, there are two definitions of the perimeter of a binary digital image. The first con- 
sists of the sum of the exposed pixel faces. To count the exposed faces we use the function 
PixelBorderLength. This takes a binary image of dimensions {3,3} , and counts the number of 
black 4-connected neighbors of the center pixel. 

Here is an example. 

Framed [test = Image [{ {0, 0, 0}, {0, 1, 1}, (1, 1, 1}}, "Bit"]] 
PixelBorderLength [test] 



2 

In this example, the center pixel has only two black neighboring pixels. 

The total perimeter can be obtained by applying this function to the {3,3} neighborhood of every 
white pixel in the image. Pixels in the interior will return a value of 0, so only perimeter pixels 
will contribute. 

To implement this idea, we first identify the positions of all the white pixels. 

positions = Position [ ImageData [jun ] , 1] ; 

Next we extract all the {3,3} neighborhoods. 

neighborhoods = 

ImageTake [ jun. Sequence @@ Transpose [ { it - 1 , tt+1}]] & /©positions; 

We can look at the first five. 
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Framed /@ neighborhoods [ [Range [5] ] ] 



Framed /@ neighborhoods [ [Range [5] ] ] 



We can also compute the pixel border length of the first five. 

PixelBorderLength /@ neighborhoods [ [Range [5] ] ] 

{ 2 , 1 , 2 , 2 , 0 } 

In a large complex image, the same neighborhood might occur many times, so we perform a tally. 

tally = Tally [neighborhoods] ; 

We can take a look at the tally, along with the pixel border length for each type of neighborhood. 
Note that 850 cases consist of all white, drawn from the interior of the foreground, with 0 border 
length. Here we just look at the first 10 elements of the tally. 
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{Framed [tt [ [ 1] ] ] , tt [ [2] ] , PixelBorderLength [tt [ [1] ] ] } & /@ tally // 
Take[tt, 10] & // TableForm 



2 


1 


2 


2 


0 


0 


0 


1 


2 


2 


The total perimeter length will be the sum of all of the border lengths for all the collected neighbor- 
hoods. We use the tallied neighborhoods, so that the pixel border length of each type of neighbor- 
hood is computed only once. 

Total [ (tt [ [2] ] PixelBorderLength [tt [[ 1] ] ] & /@ tally)] 

606 

Thus the complexity would be 

606 * 2 / 1467 / (4 Pi) // N 

19.9207 
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m Using “PolygonalLength” 

The second definition of the length of the perimeter is the sum of sides of a the polygon defined by 
the perimeter pixels considered as points in a lattice. 

To extract the perimeter, we use a new function Perimeter, defined above. This duplicates a 
builtin Mathematica function. We verify that they yield the same results. 

perimeter = Perimeter [jun] 



perimeter2 = MorphologicalPerimeter [ jun, CornerNeighbors -> False] 



ColorCombine [ { jun, perimeter, perimeter2}] 



We identify the positions of all of the pixels in the perimeter. 

Length [positions = Position [ImageData [perimeter] , 1] ] 

498 


We extract all of the 3x3 neighborhoods of pixels in the perimeter. 

neighborhoods = 

ImageTake [perimeter , Sequence @@ Transpose [ { it - 1 , tt+1}]] & /@ 
positions ; 

Then we use a new function PixelPathLength. This looks at a 3x3 binary neighborhood, identi- 
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fies the two closest white pixels not at the center, and finds the distances from their centers to that 
of the central pixel. That is the path length corresponding to that neighborhood. 

We apply this to all the perimeter pixels, and add up the results. 

Total [PixelPathLength /@ neighborhoods] // N 

547.084 

We can verify this is the same perimeter length obtained from PerimetricComplexity using Mathe- 
matica’s morphological operators. 

PerimetricComplexity [jun, Filter -> None, 

Method "PolygonalLength" ] // N 

{547.084, 1467., 16.2356} 
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